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ABSTRACT 

Context. The spectrum of IRAS 16594-4656 shows shock-excited H2 emission and collisionally excited emission lines such as [Oi], 
[C i], and [Fell]. 

Aims. The goal is to determine the location of the H2 and [Fe 11] shock emission, to determine the shock velocities, and to constrain 
the physical properties in the shock. 

/Wefhods. High resolution spectra of the 1-0 S(l), H2 2-1 S(l), [Feii], and Pa/3 emission lines were obtained with the near infrared 
spectrograph Phoenix on Gemini South. 

Results. The position- velocity diagrams of H2 1-0 S(l), H2 2-1 S(l), and [Feii] are presented. The H2 and [Fe 11] emission is spatially 
extended. The collisionally excited [O i] and [C i] optical emission lines have a similar double-peaked profile compared to the extracted 
H2 profile and appear to be produced in the same shock. They all indicate an expansion velocity of ~8 km s"' and the presence of 
a neutral, very high-density region with about 3 x 10* to 5 x 10^ cm"''. However, the [Feii] emission is single-peaked. It has a 
Gaussian FWHM of 30 km s"' and a total width of 62 km s"' at 1 % of the peak. The Pa/? profile is even wider with a Gaussian FWHM 
of 48 km s"' and a total width of 75 km s"' at 1% of the peak. 

Conclusions. The H2 emission is excited in a slow 5 to 20 km s"' shock into dense material at the edge of the lobes, caused by the 
interaction of the AGB ejecta and the post-AGB wind. The 3D representation of the H2 data shows a hollow structure with less H2 
emission in the equatorial region. The [Fe 11] emission is not present in the lobes, but originates close to the central star in fast shocks 
in the post-AGB wind or in a disk. The Pa/? emission also appears to originate close to the star. 

Key words. Line:profiles - Shock waves - Stars: AGB and post-AGB - Stars: winds, outflows - Stars: individual : IRAS 16594-4656 
- ISM: molecules - Infrared: ISM, stars 



1. Introduction 

Post- Asymptotic Giant Branch (AGB) stars represent an impor- 
tant transition phase in the evolution of low and intermediate- 
mass stars, between the AGB and the planetary nebula (PN) 
phases. During this period, the detached circumstellar envelope 
of gas and dust is expanding away from the star Meanwhile the 
star itself is increasing in temperature at nearly constant lumi- 
nosity. This phase lasts a few hundred to a few thousand years 
depending upon the star's core mass. When the temperature is 
high enough and the star photoionizes the nebula, it has entered 
the PN phase (e.g. Kwok 1993). In spite of extensive study, the 
evolution from the AGB toward the PN stage is still poorly un- 
derstood. The drastic changes observed in circumstellar structure 
and kinematics are particularly puzzling. During the late AGB or 
early post-AGB evolutionary stages, the geometry of the circum- 
stellar material changes from more or less spherically symmet- 
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ric to axially symmetric, with the result that most PNe exhibit 
axisymmetric structures, ranging from elliptical to bipolar (e.g., 
BaUck & FrankEOOlll. 

Bipolar PNe frequently possess molecular envelopes that are 
readily detectable in the near-infrared ro-vibrational lines of H2 
(Kastner et al. [T996l Kelly & Hrivnaki2005]l. The available data 
suggest that the onset of near-infrared H2 emission in PNe can 
be traced back to the proto-planetary nebula (PPN) phase but 
not back to the AGB phase of evolution (Weintraub et al. 119981 
Davis et al. |2005l l. These observations suggest that further stud- 
ies of H2 emission from PPNe may offer insight into the transi- 
tion from AGB star to PN and from spherical to axisymmetric 
mass loss. The study of transition objects showing H2 emission 
at an early stage is crucial for understanding the hydrodynamic 
processes shaping the nebulae. The H2 lines can reveal details 
about the physical conditions in the shocks associated with these 
hydrodynamic processes and thereby help constrain models of 
the interaction of the central star with the AGB remnant. 

IRAS 16594-4656 is classified as a post-AGB star for sev- 
eral reasons. In the IRAS color-color diagram it has the colors 
of a PN (Van de Steene & Pottasch il993 ), but it has not been 
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detected in the radio continuum at 3 or 6 cm. It has a large in- 
frared excess due to dust with a color temperature of 173 K. It 
displays a double-peaked spectral energy distribution, with the 
peak in the mid-infrared much brighter than the peak in the near- 
infrared (Van de Steene et al. l2000ab . It possesses a CO envelope 
with an expansion velocity of 14 kms ' (Woods et al. I2005I I. 
The chemistry of IRAS 16594-4656 appears to be carbon-rich. 
This is based on the detection of unidentified IR emission fea- 
tures at 3.3, 6.2, 7.7, 8.6, 11.3, 12.6, and 1 3.4 //m, as well as 
the 21 fim feature (Garcfa-Lario et al. 1 19991 1. all commonly as- 
sociated with a carbon-rich chemistry. The optical spectrum of 
IRAS 16594-4656 shows a spectral type B7 with significant 
reddening. The Ha emission has a P-Cygni type profile indica- 
tive of a stellar wind (Van de Steene et al. l200 0b ). In Van de 
Steene & van Hoof (120031 hereafter Paper I) we derived the to- 
tal extinction value (i.e., interstellar and circumstellar extinction 
combined) using the extinction law of Cardelli et al. ( 119891 ). We 
found Ay - 7.5 + 0.4 mag with Ry = 4. With this extinction 
value and the flux calibration from the Kurucz model, we deter- 
mined a distance of (2.2 + 0.4) L4'''^ kpc, which was in good 
agreement with the distance of 2.5 L4''^ kpc derived by Su et al. 
( 1200 11 1 (with L4 the luminosity in units of 10"* Lq). 



In the optical the nebulosity is dominated by scattered light. 
HST WFPC2 images showed that the object is a multipolar re- 
flection nebula with 3 extensions on each side (e.g., Hrivnak et 
al. II999I 1. It has been interpreted as a bipolar outflow episodi- 
cally channeled from a rotating/precessing torus. The outer neb- 
ula also shows concentric arcs (Hrivnak et al. 2001 ). High reso- 
lution mid-infrared images in the N and Q-band show a bright 
equatorial torus viewed almost edge-on and a pair of bipolar 
lobes, which show a close correspondence with the H2 HST map. 
The shape of the bipolar lobes shows that the fast outflow is still 
confined by the remnant AGB shell (Volk et al. 2006). The elon- 
gation of the inner nebula is along the symmetry axis of this 
torus at ~ 80° PA, and it seems likely that the material ejection 
is channeled into the bipolar lobes by the torus and/or some other 
collimation mechanism (Garcia-Hemandez et al. |20041 l. Ueta et 
al. (2005 2007 ) found, based on polarized-flux images, that the 
inner structure of the nebula is clearly elongated in the east- west 
direction (5'.'0 x 2'.'0 at PA 81°). Their polarized flux maps un- 
cover the bipolar cusp structure of the shell, corresponding to the 
main elongation in intensity, and a hollow shell structure which 
delineates the wall of the elongated bipolar cavities. 



In Paper I we examined the near-infrared spectrum of IRAS 
16594-4656. It shows strong H2 emission lines and some typical 
metastable shock-excited lines such as [Fe 11] 1.257 & 1.644 fim. 
We argued that the molecular hydrogen emission is mainly colli- 
sionally excited in a C-type shock. However, H2 and [Fe 11] emis- 
sion don't usually coexist. The goal of this paper is to clarify the 
location of the H2 and [Fe 11] shock emission in this object and 
to determine the shock velocities. 



High resolution near-infrared spectra of IRAS 16594-4656 
have been obtained with the near infrared high resolution spec- 
trograph Phoenix on Gemini South. The observations and data 
reduction are discussed in Section |2] The results of the analy- 
sis are presented in Section [3] and the morphology, kinematics, 
and shock properties are discussed in Section]?] and finally we 
present some conclusions in Section]5] 



2. Observations and data reduction 

2.1. High resolution infrared spectra 

The spectra were obtained in service mode with the near infrared 
high resolution spectrograph Phoenix (Hinkle et al. 12002b on 
Gemini South at 3 position angles (PA), through regions show- 
ing bright H2 emission : at 7°, 65°, and 125°. The position angle 
is measured in degrees east of north. The slit positions are shown 
on top of the HST H2 1-0 S(l) image in Fig. [T] The sHt length 
is 14" and slit width is 3 pixels or 0'.'25, giving a resolution of 
-60,000. On 5 May 2003 the integration times were 4 x 300 s 
for H2 1-0 S(l)2.12125 jum through filter K4748/4712.82 and 
4 X 200 s for [Feii] 1.64355 yum through filter H6073/6082.93 
at each position angle. H2 2-1 S(l) 2.2471 //m was observed for 
4 X 300 s at each position angle through filter K4484/4448.98 
on 1 August 2003. Next, in the same night. Pa/? was observed, 
however the weather conditions were quickly deteriorating and 
we obtained only one useful image of Pa/? before the dome had 
to be closed. 

In the archive there was another set of H2 1-0 S(l) data avail- 
able for this object taken with Phoenix under GS-2003A-Q-27 
(PI: B. J. Hrivnak; Hrivnak et al. 2006 ) at 33°, 52°, 72°, 345° PA 
through the star, and at 345° PA offset 0'.'9 west of flie star, all 
with a slit width of 0'.'35 giving a resolution of ~ 50, 000. We 
include this data set in our analysis for the sake of completeness 
and in order to reconstruct the 3-dimensional structure of the H2 
shefl. 

Basic data reduction (such as flatfielding and sky subtrac- 
tion) was done in iraf using the scripts provided by the Phoenix 
team. Next the images were rectified by fitting the position of 
the continuum peak in the spatial dimension along all pixels 
in the spectral dimension. Then the images were combined to 
generate cross- correlated, co-added images. Continuum emis- 
sion was subtracted from the data by fitting the emission of the 
source in the spectral dimension before co-adding. The wave- 
length calibration was done using telluric lines in the standard 
star spectrum using the HITRAN 2000 database (Rothman et al. 
2003 ). The resulting position-velocity diagrams are presented in 
Figs.|2]-|4] 

2.2. H2 1-0 S(1) image 

We retreived the high-resolution H2 image from the HST 
archive. IRAS 16594-4656 was observed with NICMOS 
(Malhotra et al. |2002] | on-board HST in Cycle 1 1 (General 
Observer program 9366, PI: B. Hrivnak) on 12 January 2003. 
The observations were made with NICMOS2, which provides 
a 19'.'2 X 19'.'2 field of view with 0'.'075 pixel"' scale, in con- 
junction with the F212N (H2) and F215N (H2 and Bry contin- 
uum) filters. We reduced the data using the the standard set of 
NICMOS calibration programs provided in iraf/stdas version 

3. l[^ A detailed account of the reduction procedure can be found 
in Ueta et al. (120051) . The reduced image is shown in Fig.[T]with 
the slit positions of the available Phoenix data. 

2.3. LIVES spectrum 

A high resolution echelle spectrum was obtained with UVES on 
the VLT at ESO Paranal in May and June 2000. A detailed ac- 
count of the observations, reduction procedure, and spectrum 
can be found in Reyniers (12002b . 

' STSDAS is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA 
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Fig, 1. HST/NICMOS H2 1-0 S(l) image (north is up and east to 
the left) showing the slit positions. The slit positions displayed 
are at 7°, 65°, and 125° PA with 0'.'25 width (this work; soHd 
line) and at 33°, 52°, 72°, 345° PA, and 345° PA offset a.'9 west 
of the star with (J.' 34 width (Hrivnak et al. 120061 dashed line). 

3. Results 

3.1. H2 1-0 S(1) 

The H2 1-0 S(l) position- velocity diagrams in Fig.|2]can be un- 
derstood in the context of the detailed morphology of H2 seen 
in the H2 NICMOS image. The slit positions are shown on the 
HST H2 image (Fig. [TJ. The slit position at 7° primarily goes 
through the center, which is mostly washed out in the NICMOS 
H2 image due to the bright central star The slit position at 65° is 
almost along the major axis and going through the bright spot at 
the tip of the eastern lobe. The slit position at 125° goes through 
the center at an intermediate angle and captures the bright north- 
western edge of the lobe. 

Reyniers (2002) observed a large spread in radial veloc- 
ity of the optical lines (from -10 to -40kms"') with a 
small dependence upon excitation potential. The high excita- 
tion emission lines seem to have mean radial velocities around 
-30kms"', while the absorption lines have velocities more 
around -23kms ' in the mean. The extracted H2 profiles are 
double-peaked. The minimum flux in the H2 profiles is on aver- 
age at -29.2 km s"'. Hence the H2 velocity seems to agree on 
average with the velocity of the emission lines in the optical 
spectrum. We adopt -29.2 km s"' as reference velocity. 

The slit position at 7° PA (Fig. |2] top left) shows a toroidal 
structure of 2'.'4 in diameter in the north-south direction at 10% 
of the peak. The diameter of the shell extends to about 4" at 
1% {~ 5cr) of the peak emission. There is little emisson around 
the central star, reinforcing the toroidal appearance. It may cor- 
respond to the toroidal equatorial density enhancement (EDE) 
as seen in the mid-infrared images (Volk et al. 2006). The bulk 
of the near/approaching side of the shell moves at -38.2 km s"' 
and the far/receding side at -21.2 km s"'. Hence, with respect to 
the reference velocity of -29.2 km s ', the bulk of H2 emission 



appears to expand at about 8kms . This seems low compared 
to the CO (J=2-l) outflow velocity of 14km s ' (Woods et al. 
120051 ). However, the total velocity width, including the wings, 
is as large as 50.8 (±2.2) km s"' at the 1% level. The velocity 
spread in the wing at the receding side is much larger than at the 
approaching side. 

The slit position at 65° PA (Fig. [2] top right) samples the 
lobes close to the major axis. In this spectral image the object 
is most elongated: about 5.5" at the 1% (2cr) level of the peak. 
The bulk of the emission is detected within 2" of the central 
star, most of it towards the northeast, which corresponds to the 
main emission region in the HST H2 image (Fig. [T]). The lack 
of emission at the reference velocity close to the star points to a 
hollow nature of the shell. At this slit position, the main lobe is 
separated into two emission regions that are moving away from 
the central region. The approaching peak is located almost at 
the center at -39.5 km s ', while the receding shell is slightly 
off (0'.'3) to the southwest and moves at -20.6 km s"'. The third 
emission peak represents the northeastern "tip" at V.'4 from the 
centre and moves at -26.7 km s ', which is close to the refer- 
ence velocity. The approaching peak is strongest and appears to 
be elongated towards the southwest. The presence of the blob 
at the northeastern tip seems to make the receding side appear 
to be elongated towards the northeast, giving a point-symmetric 
appearance to the whole shell structure. 

At 125° PA (Fig. |2] bottom second left) the slit still passes 
through the EDE, but also through the region of bright H2 emis- 
sion in the western lobe. The shell structure appears to be hollow, 
which is consistent with what is seen in the data at other slit po- 
sitions. The position-velocity diagram suggests that most of the 
bright emission at the northwestern side of the lobe comes from 
the approaching side of the lobe (at 0'.'7 to the northwest with a 
velocity of -36.8 kms"') 

Data at other slit positions have been discussed by Hrivnak 
et al. (2006). The overall trend of the emission structure is sim- 
ilar for the two data sets at similar position angles, except that 
the Hrivnak data set covers emission from the clumps located at 
~ 5" from the centre. The clumps are almost at the reference ve- 
locity. However, especially at the 52° PA, there is a slight veloc- 
ity difference between these clumps: the clump at the northeast 
is a bit blueshifted and the one at the southwest is a bit redshifted 
compared to the reference velocity. The slit positions at 33° and 
125° PA, and as well as at 72° and 345° PA are almost perpendic- 
ular to each other. The data confirm the hollow structure which 
is round and elongated, but the H2 intensity is not uniformly dis- 
tributed along the walls. 

In general the position- velocity diagrams suggest that (1) the 
shell is hollow, (2) its axis is nearly aligned with the plane of 
the sky, (3) the approaching side of the shell shows stronger 
emission and is elongated towards the southwest, while the posi- 
tion of the emission peak is slightly northeast of the centre, and 
(4) the receding side of the shell is weaker in emission strength 
and is elongated towards the northeast, while the position of the 
emission peak is slightly southwest of the centre. Thus the emis- 
sion structure appears to be almost point-symmetric in both the 
spatial and velocity dimensions. 

3.2. H22-I S(1) emission 

Fig. [3] shows the H2 2-1 S(l) position-velocity diagrams (left 
column) and H2 2-1 S(l) contours overlaid on the H2 1-0 S(l) 
position-velocity diagrams (right column at 7°, 65°, and 125° 
PA, from top to bottom, respectively). The H2 2-1 S(l) and 1-0 
S(l) emission seem to be extended in a similar manner, except 
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Fig. 5. Velocity profiles of the extracted spectra of H2 1-0 S(l) 
at 7° PA, [Fe 11] at 7° PA, [O i], and [C i]. Each line is continuum 
subtracted and normalized to the peak intensity. 



Fig. 6. Velocity profiles of the extracted spectra of 3 Paschen 
lines and [Fen]. Each line is continuum subtracted and normal- 
ized to the peak intensity. 



that (1) H2 2-1 S(l) seems to originate a bit closer to the cen- 
tral star than H2 1-0 S(l) and (2) H2 2-1 S(l) does not show an 
elongation towards the northeast (in the receding lobe) and to the 
soutwest (in the approaching lobe) at 65° PA. Also, H2 2-1 S(l) 
does not seem to show the large velocity dispersion towards the 
red as in the H2 1-0 S(l) images. The lowest contour (at 1% of 
the peak) displays elongations into higher/lower velocities near 
the centre, but this is probably residual continuum emission that 
has not been completely removed. 

3.3. [Feu] 1.645jdm emission 

Fig. |4] shows the [Feii] position-velocity diagrams (left column) 
and the [Feii] contours overlaid on the H2 1-0 S(l) position- 
velocity diagrams (right column) at 7°, 65°, and 125° PA, from 
top to bottom, respectively. From the contour-overlaid images it 
is clear that the H2 and [Fe 11] emission originates from different 
regions. 

The seeing of the night was between 0'.'7 and 0'.'9 in the 
R-band. Using the stellar continuum we measure a seeing of 
a'8, 0'.'7, and 0'.'6 on average at 7°, 65°, and 125° respectively. 
The Gaussian FWHM of the [Feii] profile is a'98, 1'.'05, and 
0'.'94 at 7°, 65°, and 125° respectively. Hence we conclude that 
the [Fell] emission region is extended. At an assumed lumi- 
nosity of 10^ L0 and distance of 2.2 kpc (Paper!) the decon- 
volved Gaussian FWHM diameter of the emitting region would 
be 0.004 pc ± 0.001 pc (or 80 + 20AU). Assuming a stellar tem- 
perature of 14,000 K (Revniers l2002i l. this is about 1800 stellar 
radii. 

The extracted profile at 7° PA is shown in Fig.|5] The [Fe n] 
emission peaks at the reference velocity of -29.2 km s"' and the 
Gaussian FWHM of the profiles are 31.8kms 33.1 kms ', 
and 33.1 kms ' at 7°, 65°, and 125° PA, respectively. The to- 
tal width of the profile at 1 % of the peak is 62 km s ' . The [Fe 11] 
emission appears to come from the inner region, close to the star 
at all position angles and shows a larger velocity dispersion than 
the H2 emission. 



3.4. Pa/3 and otiner Pascinen lines 

One useful image of Pa/? was obtained at 7° PA before the dome 
had to be closed due to bad weather The Pa/3 emission is cen- 
tered on the star. The seeing was at least r.'5, much worse than 
for the [Fe 11] observations, and the spatial profile remained un- 
resolved. 

In Fig.|6]we compare Pa^S with the [Feii] profile at 7° PA. 
The peak of the Pay6 profile is at -21 km s while the peak of 
[Fell] is at -29.8 kms The velocity of the peak of the Pa^S 
profile is closer to the average velocity of the absorption lines 
in the optical spectrum, while the velocity of the peak of the 
[Fe 11] emission is closer to the velocity of the the optical emis- 
sion lines. Pa^ is faint, but the normalized profile is consider- 
ably broader than the [Feii] emission profiles. The Gaussian 
FWHM is 48.2 kms"' and the total profile width at 1% of the 
peak is ~75kms"'. The peak of the profile appears to be asym- 
metric. In Fig. [6] we compare the Pa/3 profile with two other 
Paschen lines (Pad at 1004.9373 nm and Pa20 at 839.2396 nm) 
from the UVES spectrum. The absorption profile of Pa 3-20 has 
very broad wings and some emission. This emission becomes 
stronger for lower members of the Paschen series. 

3.5. Optical [Oi] and [Ci] emission lines 

The only forbidden emission lines in the optical, UVES spec- 
trum (Reyniers|2002ll are the [Oi] lines at 5577.3387, 6300.304, 
6363.777A and the [Ci] lines at 8727.12, 9824.13, and 
9850.26 A. We have plotted the strongest line of each element 
in Fig.|5] We notice that the optical [O i] and [C i] emission lines 
have a profile very similar to the H2 emission line. Because the 
UVES spectra were obtained with the image sheer, we have no 
spatial information. However based on their similar profiles, they 
seem to originate in the shock which produces the H2 emission. 
In PNe, collisionally excited C*^ has also been observed to co- 
exist mainly with H2, making it an invaluable probe of H2 emit- 
ting regions (Liu et al. 1995). In old extended PNe, such as the 
Ring nebula, the [O i] 6300 A line is also found to largely follow 
the near-infrared H2 emission (Liu & Barlow rT996b . 
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The [C i] A 9824, 9850 A lines decay from the same upper 
level, so their intensity ratio depends only on the relative transi- 
tion probabilities, and is purely determined by atomic parame- 
ters. Similarly, the [Oi] /1 6300, 6363 A line ratio is also purely 
determined by atomic parameters. In the case of collisional exci- 
tation, the observed [C i] line ratio (79824 + 79850) / 78727 and 
the [O i] ratio (76300 + 76363) / 75577 provide each a direct di- 
agnostic of the electron temperature, assuming the electron den- 
sity is known. Using the extinction derived in Paper I (Ry - 4.2, 
E(B-V) = 1.78), we corrected the lines for extinction and cal- 
culated both ratios. For the [Ci] line ratio we obtain a value of 
1.5 ±0.5. For the [Oi] ratio we obtain a value of 6.0± 1.0. The 
UVES slit width was 0'.'7, orientated east-west, and, because of 
the use of the slicer, no atmospheric dispersion correction could 
be applied. The shift between the 5577 A and the 6363 A lines is 
0'.'2. Because the object is extended, this probably has little effect 
on the ratio. We checked whether the [Ci] lines are affected by 
telluric absorption using the transition probabilities calculated 
by Nussbaumer & Rusca ( fT979] ). The 79824 / 79850 theoretical 
line ratio is 0.338 and our value is within 3% of this theoretical 
value. Hence in our spectra these lines are not affected by telluric 
absorption 

In order to calculate the electron temperature from the [O i] 
or [Ci] line ratios, we need to assume an electron density. 
Usually, the [Ni] doublet ratio A 5200.42/5197.95 A is one of 
the very few density-sensitive diagnostics observable from neu- 
tral species in the optical region. It is a good tracer of the electron 
density in regions where [C i] and [O i] are observed. However 
the [N i] lines are not detected. This is not due to ionization of ni- 
trogen. The ionization potential of N" (14.534 eV) is a bit larger 
than that of H° (13.62 eV) and O" (13.60 eV), and no [On] 
or [Nil] lines are observed in the optical spectrum. Nitrogen 
doesn't seem to be underabundant (Revniers l20021 l either. Hence 
the absence of the [N i] lines indicates that the [N i] lines are col- 
lisionally de-excited in the high-density region where the [C i] 
and [Oi] lines originate. This is due to the relatively low crit- 
ical densities of the [Ni] lines: about 4830 and 1160 cm"-' at 
Te = 10, 000 K for the upp er lev els of the i 5198 an d 5200 A 
lines, respectively (Zeippen i 19821 Berrington & Burke [T98Tb . 

For carbon, the critical densities of the 2p^ 'D2 level from 
which the /1 9824 A and 9850 A transitions originate are 4.8 x 
lO'*, 2.2 X lO'*, and 1.6 x 10^ cnrffor = 1000, 5000, and 
10,000K, respectively (Liu etal.|T995ll and 1.1 x lO^cm^^for 
the 'So level from which the /1 8727 A line originates (Mendoza 
119831 Berrington & Burke [T98T ). Hence the low value derived 
from the [Ci] ratio indicates that the 'D2 level population is 
quenched by electron impacts, and that the electron density in 
the shocked region would be higher than 10^ cm"^. The crit- 
ical density for the 'D2 level from which the /1 6300 A and 
6363 A lines originate is 1.9 x lO^'cm"^ at 1 0,000 K and is 
1.1 X 10*^ cm"^ for the 'S q level from which t he /1 5 577 A line 
originates (Mendoza [T983l Berrington & Burke 19811. The [Oi] 
line ratio is insensitive to electron densities below approximately 
lO*" cm-^ at 10,000 K. 

As we have no constraint on the electron density, we have 
used the photo-ionization code Cloudy (version 07.10.15) last 
described by Ferland ( 1998 ) to calculate a grid of [O i] and [C i] 
line ratios in function of and tif.. The detection of strong col- 
lisionally excited [C i] and [O i] line emission, but no [0 11] line 
emission, requires that there must be a region which is predomi- 
nantly neutral, but hot enough for effective excitation of the lev- 
els in [Oi], which implies 6000 ^ re(K) < 10,000. The elec- 
tron density must be higher than 10"* cm"^ for the [N i] lines to be 



collisionally de-excited. We considered 3.7 :$ log(re/K) ^ 4.0 
and 4.0 :^ log(ne/cm"-') ^ 8.0, using steps of 0.1 dex and 0.25 
dex respectively. The [O i] ratio corresponds to electron densi- 
ties of 3 X 10^ to 5 X 10^ cm"-' for electron temperatures of 
10,000 K to 6,000 K, respectively. At these electron tempera- 
tures the electron densities derived from the [C i] line ratio are 
10^ - 3 X 10^ cm-^ 

In summary: There is no indication of ionized gas in the opti- 
cal spectrum of IRAS 16594-4656. The forbidden [O i] and [C i] 
emission are collisionally excited and originate in a very high 
density region (3 x 10^ :$ :$ 5 x 10^ cm"-'), in the same 
shock which produces the H2 emission. 

4. Discussion 

4.1. Kinematics and morphology of the H2 1-0 S(1) emission. 

In order to investigate the 3-dimensional structure of the lobes, 
we constructed a data cube using all existing data in H2 1-0 S(l). 
First we rectified the H2 1-0 S(l) data at each position angle 
to have a pixel scale of 0.085" pix"' in the spatial dimension 
and of 1.263 km s"' in the velocity dimension. Then we aligned 
the position-velocity frames using the location of the contin- 
uum peak in the spatial dimension and the reference velocity 
in the velocity dimension. Next we spline-interpolated the sur- 
face brightness as a function of position angle in each pixel in 
the position velocity frame based on the original seven frames. 
Finally we projected the interpolated 3-dimensional data onto a 
3-dimensional cartesian grid. 

Fig.|7] shows the velocity channel maps that are recovered 
from the reconstructed data cube in H2 1-0 S(l). Each channel 
has a width of 5.05 km s"'. Maps near the reference velocity in- 
dicate that the H2 emission is strong in both lobes, but stronger 
in the western lobe than in the eastern lobe, which is also seen in 
the total intensity map in Fig.[T] They also indicate that there is 
in fact less emission along the waist region of the shell. This can 
also be seen from the position velocity maps at 7° and 345°, but 
the velocity channel maps depict it more eloquently. Thus, the 
overall 3-dimensional structure of the shell is such that the ma- 
jority of the H2 emission is distributed as hollow bipolar lobes, 
one at the approaching side and one at the receding side of the 
shell. The approaching side of the shell is elongated towards the 
west while the receding side of the shell is elongated towards the 
east. 

Assuming a homologous expansion, i.e., an expansion with 
V oc r, the velocity axis in the recovered data cube corresponds 
to a spatial axis along the line of sight. This allows us to vi- 
sualize the 3-dimensional structure of the spatial distribution of 
H2. Fig.[8]shows the 3-dimensional representation of the central 
shell of IRAS 16594-4656, for which the surface contour has 
been constructed based on the reconstructed data cube at 30% 
of the peak emission. Fig. [8] shows the bottom right view of the 
surface contour, illuminated obliquely from the southwest side 
by a fiducial light source. The 3-dimensional volume shows the 
central cavity and the hollowness of the shell along the polar axis 
of the bipolar lobes, the western and eastern elongations in the 
approaching and receding sides of the shell, and the holes on the 
surface at the top and bottom of the equatorial regions. 

4.2. Shock properties 

In Paper I we concluded that the H2 emission in IRAS 
16594-4656 is mainly collisionally excited and that (1) the 
line ratios indicate that the H2 emission originates in a 20 to 
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Fig. 8. 3-dimensional rendition of the shell structure of IRAS 
16594-4656. The surface represents an isointensity contour at 
roughly 30% of the peak intensity. 



30km s ' C-type shock in 10^ cm"^ material (Le Bourlot et al. 
I2OO2I 1. (2) we postulated that the C-shocks occur where the fast 
wind is funneled through the EDE, which may harbor a magnetic 
field, but we noted that the H2 emission could also be present in 
the lobes, excited by shocks caused by a molecular outflow im- 
pinging on the AGB envelope, and (3) we noted that in principle 
H2 and Fe^ cannot coexist in the same shocked region in sub- 
stantial quantities. 

In previous sections we showed that the H2 emission is in- 
deed present at the edge of the EDE, but is more intense at the 
edge of the lobes. The H2 emission appears to be excited in 
shocks between the AGB ejecta and the post- AGB wind. 

In this paper we found that the [O i] and [C i] emission lines 
are very similar to the H2 profiles in shape and velocity. They 
may originate in the shock close to where the H2 emission orig- 
inates. Their line ratios indicate a very high density of 3 x lO*" 
to 5 X 10^ cm"^ in the shock. At this high density the rotational 
temperature of 1440+ 80 K derived in Paper I corresponds to 
a 10 to 20 kms"' C-type shock into lO*" to 10^ cm"^ material 
(Le Bourlot et al. 2002). However at such high densities, the 
ionization of the metals may suffice to produce a J -type shock 
(Smith [r994l) . If this were the case, the rotational temperature 
determined in Paper I and the high density would also be consis- 
tent with a J-type shock of 5 to 10 km s"' into 10^ to 10^ cm"-' 
material (Wilgenbus et al. 2000). At these low velocities and 
high densities the J-type shock will not dissociate H2, and the 
H2 emission will remain strong. 

[Fell] 1.64 fim is a cooling line of gas that is excited pre- 
dominantly by shocks. It is a very good tracer of gas shocked 
by energetic events in the absence of photoionization, which is 
the case for IRAS 16594-4656. [Feii] emission cannot be ex- 
cited in the presence of C-shocks alone, as C-type shocks do not 
produce emission from ionized species such as Fe^. [Fe 11] emis- 
sion is well modelled by fast dissociative J-type shocks. In slow 
non-dissociative J-shocks, the [Fe 11] emission drops rapidly with 
the shock velocity (Gredel jl994 ). Hence [Feii] emission arises 



from fast dissociating and ionizing shocks, while the H2 emis- 
sion arises from locations where the shock velocities are low 
enough, not to cause destruction of H2. We see indeed that the 
H2 and [Fe 11] emission line profiles and spatial extent are suf- 
ficiently different to confirm that they mainly come from differ- 
ent regions. The [Feii] emission appears to originate in a region 
close to the central star within the hollow structure outlined by 
the H2 emission. The [Fe 11] emission is not located along the 
side walls of the hourglass, as was considered as one of the pos- 
sibilities in Paper I and is the case in Hubble 12. Neither is it 
present as blobs in discrete shock fronts in the lobes as is the 
case for M 1-92 (Davis et al. 2005). The [Fen] emission may be 
due to shocks in the stellar wind or originate in a circumstellar 
or circumbinary disk. 

5. Conclusions 

The H2 emission is excited in slow 5 to 20 kms"' shocks into 
dense material at the edge of the lobes, caused by the interac- 
tion of the AGB ejecta and the post-AGB wind. The 3D repre- 
sentation of the H2 emission shows a hollow structure. There is 
less H2 emission in the equatorial region. The collisionally ex- 
cited [Oi] and [Ci] optical emission lines have a similar profile 
compared to the extracted H2 profile and appear to be produced 
in the same shock. They all indicate an expansion velocity of 
~ 8 km s ' and the presence of a neutral, very high density re- 
gion of about 3 X 10^ to 5 x 10^ cm"-'. The [Feii] emission is 
not present in the lobes, but originates close to the central star. It 
originates in fast shocks in the post-AGB wind, or in a circum- 
stellar or circumbinary disk. The Pa/? emission also appears to 
originate close to the star. The total width of the [Fe 11] and Pa/3 
emission line profiles is 62 and 75 km s"' respectively, at 1% of 
the peak. 
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Fig. 2. The position-velocity diagrams of the spectra at H2 1-0 S(l). The slit position angles and corresponding spatial directions are 
indicated at the top and bottom comers. Contours are, from the top to bottom, 90%, 70%, 50%, 30%, 10%, 5%, and 1% of the peak. 
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Fig. 3. The position-velocity diagrams of the spectra at H2 2-1 S(l) (left) and H2 2-1 S(l) contours overlaid on the H2 1-0 S(l) 
images (right). The slit position angles and corresponding spatial directions are indicated at the top and bottom comers. Contours 
are, from the top to bottom, 90%, 70%, 50%, 30%, and 10% of the peak. When conturs are overlaid, the H2 1-0 S(l) images are 
smoothed to have matching spatial and spectral resolution. 
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Fig. 4. [Fell] color images with [Feii] contours overplotted (left) and [Fen] contours overplotted on the H2 1-0 S(l) images (right) 
at 7°, 65°, and 125° PA from top to bottom, respectively. 
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Fig. 7. The velocity channel maps of the reconstructed H2 1-0 S(l) data cube centered around the reference velocity of 
Each channel has a width of 5.05 km s ' and its central velocity is indicated at the bottom left corner in km s ' . 



